Generally covariant model of a scalar field with high frequency dispersion 
and the cosmological horizon problem 
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Short distance structure of spacetime may show up in 
the form of high frequency dispersion. Although such dis- 
persion is not locally Lorentz invariant, we show in a scalar 
field model how it can nevertheless be incorporated into a 
generally covariant metric theory of gravity provided the lo- 
cally preferred frame is dynamical. We evaluate the resulting 
energy-momentum tensor and compute its expectation value 
for a quantum field in a thermal state. The equation of state 
differs at high temperatures from the usual one, but not by 
enough to impact the problems of a hot big bang cosmology. 
We show that a superluminal dispersion relation can solve 
the horizon problem via superluminal equilibration, however 
it cannot do so while remaining outside the Planck regime 
unless the dispersion relation is artificially chosen to have a 
rather steep dependence on wavevector. 



I. INTRODUCTION 

Recent studies of the role of trans-Planckian physics 
in the Hawking effect and in inflationary cosmology have 
exploited scalar fields with high frequency dispersion as 
a model of how short distance physics might affect the 
behavior of quantum fields in such settings. These mod- 
els involve a preferred frame in which the distinction be- 
tween high and low frequency is made. If such a matter 
field is to be coupled to the spacetime metric in a gener- 
ally covariant theory of gravity, the preferred frame must 
be treated as a dynamical quantity rather than as a fixed 
background structure. One way of doing this was formu- 
lated by us in a previous paper |l| , in which the preferred 
frame is determined by a dynamical unit timelike vector 
field w°. Using this formulation we obtained an expres- 
sion for the stress tensor of a scalar field with dispersion. 

In the present paper we use this stress tensor to eval- 
uate the equation of state in fiat spacetime for a thermal 
state of the field. This equation of state is then used in 
a fiuid description of the matter field in a cosmological 
model, and the implications for a hot big bang cosmology 
are examined. We also investigate under what circum- 
stances superluminal equilibration associated with high 



frequency dispersion can solve the cosmological horizon 
problem. 

In considering only the flat space thermal state of the 
matter field we are adopting an adiabatic approxima- 
tion which precludes effects related to out of equilibrium 
phenomena including particle creation in the dynamical 
background of a cosmological metric such as were con- 
sidered recently for fields with high frequency dispersion 
in the context of inflation |g,|[ . This is a good approxi- 
mation for frequencies larger than the expansion rate H . 
It can therefore be used to study the effects of dispersion 
for frequencies of order fcp provided H ^ kg. 



II. MODEL FIELD THEORY 

Various possibilities exist for the kinetic terms in the 
action for m". For the purposes of illustration in this 
paper we choose the "minimal theory" of Ref . M , 



(-aiR - h F^'^'Fab + Kgabu'^u'' - 1)) , (1) 



where 



Fab := 2V[QMf,]. 



(2) 
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The field A is just a Lagrange multiplier whose variation 
enforces the constraint that u° be a unit vector. 

For the matter content we are interested in a scalar 
field with high frequency dispersion oj'^ = \k\'^[l + 
.g(|A:|/fco)], where g is a function that vanishes at zero, 
and ko is a constant with the dimensions of inverse length 
which sets the scale for deviations from Lorentz invari- 
ance. It has been suggested that such a modified disper- 
sion relation might arise in loop quantum gravity P,p|, 
or in string theory or other approaches to quantum grav- 
ity, or more generally from an unspecified modification of 
the short distance structure of spacetime (see for exam- 
ple gj^). Possible observational consequences have been 
the subject of recent study (see for example PP-[lO|| and 
references therein), and the role of such dispersion in the 
Hawking process [Q and in the generation of infiationary 
primordial density fluctuations MM have been examined. 

Absent a reliable theory of such modifications, it makes 
sense simply to expand in fc/feg. We consider here the 
lowest order modification for a scalar field that is invari- 
ant under rotations and analytic in k, which is given by 



u;ikf = \k\'~\k\ykl 



(3) 



where fcg is a dimensionful parameter. With fcg > this 
yields a subluminal group velocity, and with fcg < it is 
super luminal. The pathology at |A:| = ko in the sublu- 
minal case is irrelevant since (||) is only regarded as the 
first terms in an expansion. When we evaluate the equa- 
tion of state, we shall impose a cutoff that avoids this 
pathology. 

The dispersion relation (ph can be produced by adding 
a term to the action with four spatial derivatives. Previ- 
ously this has been done in 1-1-1 dimensional models (for 
a review see [[7| ) , and recently such models have been gen- 
eralized to field theory in a background 3-1-1 dimensional 
Robertson- Walker spacetime [g,Q] . Here we extend these 
models to a general 3-1-1 dimensional spacetime. This can 
be accomplished in a generally covariant manner, consis- 
tently with spatial rotation invariance, with the action 
S^ ~ J d^x ^—g C^p, where 

i^^ = \ (vVVa'/' + K\D^^f) . (4) 

Here D^ is the covariant spatial Laplacian, i.e., 

Z?V = -D^Da^ = -g"=Va(gc'V6(^), (5) 

where Da is the spatial covariant derivative operator |l^] 
and Qab is the (positive definite) spatial metric orthogonal 
to the dynamical unit vector m", 



Qab ■= -gab + UaUb. 



(6) 



The equation of motion for the metric takes the Ein- 
stein form 



(jab 



inG{T^:^+Tifi, 



(7) 



where G — l/Wirai, and with 

Tat^ = -^bi{FaraFb^ - ^F^gab) + 2XUaUb, (8) 



^ah = ^a^P^b-i^ - C^ gab 



V'^'iqabD^ipDrnip) 



(9) 



(The constraint equation gabu°'u^ = 1 has been used to 
drop the contribution to (0) that would have come from 
the variation of ^—g in the constraint term of the action 

(DO 

The equation of motion for the field u" takes the form 



V'Fba 



2^K 



15S^ 



2hi\ 2 Su'' 

with 



(10) 



Contracting (|T^) with w" we obtain an expression for A 
in terms of the other fields: 



2 d-u" 



where, from (|l l|) , 



w"^ = Ak^^u'^v'D^^^^^^^a]qb 



(12) 



(13) 



III. THERMAL STATE IN A 
ROBERTSON- WALKER COSMOLOGY 

Let us now specialize to a Robertson- Walker (RW) 
spacetime, in the semiclassical framework where the met- 
ric and w° are treated as classical fields and the scalar 
field is a quantum field (which therefore has a well- 
defined thermal equilibrium state). In the field equations 
for the metric and for u° we take the expectation value 
of the (/?-terms. Assuming u°' shares the RW symmetry 
it must be the cosmological rest frame. The tensor Fab 
then vanishes, so A is just determined by the matter term 
in (|l2). (If there are further terms involving u° in the 
action then there are additional contributions to A.) 

With RW symmetry we have V aUb = Hqab, where 
H = a/ a is the usual Hubble "constant". Iterating this 
identity we find VaQmn = H{qamUn + qanUm), from 
which it follows that the contraction (|l^) is given by 



«"^=6fco-'^^^V, 



(14) 



(11) 



where ip = u"^'S7rnf- 

Suppose now that the scalar field is well approximated 
by an adiabatically evolving thermal state. This would be 
the case if (i) there are interactions that produce an equi- 
libration rate which is large compared to the expansion 
rate, and (ii) the thermal frequency is also large com- 
pared to the expansion rate. In this case the expectation 
value of (the Hermitian part of) the operator ip D'^ip van- 
ishes, since the operator is odd under time reversal while 
a thermal state is invariant. Using (|lj) and (|lj) this im- 
plies that, in the minimal model, (A) — 0. (If other terms 
are included in the action for m° then (A) would not van- 
ish, although it would still not receive contributions from 
this matter field.) Therefore in this model {T^^ ) = 0, so 
the only contribution to the cosmological stress tensor 
comes from the scalar field. 



A. Thermal equation of state 

Consider now the expectation value {T^ ) of the stress 
tensor (0) — or more precisely of its Hermitian part — in a 
thermal state. Note first that time reversal invariance of 



the thermal state requires an even number of time deriva- 
tives of (/?, and spatial isotropy requires an even number 
of spatial derivatives oi tp, in order for the thermal ex- 
pectation value not to vanish. Thus terms with an odd 
number of derivatives of ip do not contribute. Let us call 
this the "odd derivative rule" . We can use this rule to see 
that the expectation value (Cip) vanishes. Integrating by 
parts, C^ can be expressed as a term that vanishes since 
(fi satisfies its equation of motion, plus the total derivative 
of an expression involving only terms with an odd num- 
ber of derivatives of tp. The total derivative can be taken 
out of the expectation value, hence that term vanishes. 

The part of @ multipled by k^^ has three terms. 
The expectation value of the first term has a single time 
derivative, hence vanishes by time reversal symmetry. 
The expectation value of the third term is the gradient 
of an expression with three spatial derivatives which van- 
ishes by the odd derivative rule. The second term can 
be integrated by parts, and the resulting total derivative 
piece has vanishing expectation value by the odd deriva- 
tive rule, which leaves only 2q(^a"^{D'^(fVm'^b)^)- In ^ 
thermal state this must have the form AuaUb + Bqab, 
where w" is the rest frame defined by the thermal bath. 
To find A we contract with u'^u^, which yields ^4 = due 
to the factor ga™- To find B we contract with g"*" and 
divide by 3, hence this term contributes —^{D'^ipD^Lp). 
The expectation value of the stress tensor thus takes the 
form 



(Tab) = pUaUb + Pqab, 

with energy density p and pressure P given by 



(15) 



(16) 
(17) 



To evaluate the density and pressure we expand the 
field in Fourier components using the dispersion relation 
(^ and sum the contributions from the modes, weight- 
ing each by the thermal expectation value of the number 
operator e^'^^'^l^ — 1, which yields 



P = 



d^k 



uj{kf 



large fc, so in the high temperature limit only wave vec- 
tors \k\ ^ \ko\ are relevant. For such wave vectors we 
have uj{k)^ ~ l^l^'/l^ol' hence the energy density scales 
as p ex T^^'^k^ , and P ~ |p. This should be taken only 
as a qualitative indication of what might occur, since 
one would expect that for |fc| > fco further terms in the 
fc-expansion of the dispersion relation become important. 

In the subluminal case k'^ > ^ we must impose a cut- 
off so as not to enter the unphysical region where the 
dispersion relation would yield aj(fc)^ < 0. We choose 
to impose the cutoff at |A:| = /co/\/2, where uiik) at- 
tains its maximum value. This is similar to a lattice 
cutoff for which, in one dimension, the dispersion rela- 
tion is cj(fc) = (2/^) sinfc(5/2, and for which |fc| < tt/J 
exhausts all the independent modes. With this cutoff in 
place, Lo{k)/T ^ 1 for all modes in the high temperature 
limit, hence the expectation value of the number opera- 
tor tends toward T/uj{k). Thus the energy density scales 
as /9 ex Tfcp in this limit. As for the equation of state, the 
ratio of the integrals yields P/p ~ 0.174. 

The interpolation between the low and high tempera- 
ture limits can be determined by a numerical calculation 
of the integrals. In Fig. |l| we plot P/p vs. log;^Q(T/fco). 
It is seen that the equation of state smoothly connects 
the low and high temperature limits. For the subluminal 
case most of the interpolation takes place over the range 
of temperatures 10^^-^ < T/fco < lO^''-^, while in the 
superluminal case most of the interpolation takes place 
over the somewhat higher range 10^^ < T/ko < 10. 



27r3u;(fc) e'^C^-)/^ - 1 ' 



(18) 



P 



(fk LOJkf - |fc|Vfcg 

27r3t^(fc) e"(^)/'^ - 1 



(19) 



For low temperatures T <^ ko, only modes with k <^ k^ 
contribute significantly, hence we recover the standard 
result for massless radiation, P = ^p, and the energy 
density scales as p ex T''. 

Before looking at the exact temperature dependence, 
let us consider the high temperature limit. The nature 
of this limit depends on whether the dispersion is sub- or 
super-luminal. In the superluminal case k^ < 0, we can 
sensibly use the dispersion relation (ph out to arbitrarily 



FIG. 1. Thermal expectation values for P/p vs. log(r/fco) 
in the presence of high frequency dispersion. The dashed and 
solid curves are for the sub- and superluminal cases respec- 
tively. 



B. Cosmological implications 

The modified equation of state derived here can be 
used in a cosmological model where thermal matter with 
non-Lorentz-invariant dispersion acts as a source of grav- 
ity. The Bianchi identity implies the energy-momentum 



tensor is divergence-free, which for a perfect fluid implies 
p = —3H{p+p). Together with the equation of state this 
determines the evolution of the temperature. In the su- 
perluminal case local energy-momentum conservation is 
also independently implied by the matter field equations 
as usual, however in the subluminal case the dynamics is 
not self-contained because of the high wave vector cutoff. 
If the cutoff is imposed at a fixed proper wave vector, new 
modes are added as the universe expands, and the field 
dynamics does not specify into which state these modes 
are born. The assumptions of a thermal stress tensor and 
of energy- momentum conservation require that, whatever 
their birth state, they rapidly equilibrate with the rest of 
the modes. This is reasonable under the assumption that 
the system is in the adiabatic regime. 

It is interesting to ask whether high frequency disper- 
sion could have any impact on the cosmological problems 
that led to the invention of the inflationary scenario. In 
particular, could it provide a solution of these problems 
not requiring inflation? The modiflcation of the equa- 
tion of state caused by the dispersion would affect the 
quantitative details of the evolution of the scale factor 
at times when the typical wavevectors are of order kg or 
greater, which is presumably only in the very early uni- 
verse. However, as we have seen, the equation of state 
only changes from P = p/3 to P = (2/3)/9 in the super- 
luminal case and to P = 0.17p in the subluminal case. 
Neither are significant enough to qualitatively change the 
dynamics or horizon size. 



C. Superluminal equilibration and the horizon 
problem 

In the superluminal case, the fact that influences travel 
faster than light at wavevectors k ^ k^ opens up the 
possibility of solving the horizon problem via superlumi- 
nal equilibration. The coordinate distance covered by 
a wavepacket with proper group velocity Vg is Aa; = 
/ Vgdt/a. For a dispersion relation that goes as w ~ fc" at 
large wave vectors the group velocity goes as Vg ~ fc"^^. 
It is easy to show that the typical wavevector at the peak 
of a thermal distribution scales as a~^ (as long as the dis- 
persion relation is homogeneous), hence the typical group 
velocity scales as a^~", so we have Aa; ~ / dt/a". If this 
diverges then the horizon problem is solved assuming the 
framework of this model. 

With the above dispersion relation, the equation of 
state is P = (n/3)p, which according to the Einstein 
equation yields the evolution of the scale factor a{t) ex 
i2/("+3). Thus Ax - fdtt-^''/^''+^\ which diverges at 
the lower limit for any 71 > 3. 

Unfortunately it does not make much sense to view 
this as a solution to the horizon problem, because most 
of the Aa; is traversed during a regime in which the typ- 
ical wavevector and energy density are so much larger 
than the Planck scale that we have no reason to trust 



the semiclassical model at all. In particular, we have 
checked that in order for Aa; to surpass the horizon size, 
the evolution must be extrapolated all the way back to a 
time at which the typical wavevector is roughly fcg times 
[10^'^(A;o/fcp;ancfc)^]^^^"~^^- Presuming that fco is within 
a few orders of magnitude of the Planck scale, this typi- 
cal wave vector exceeds the Planck scale unless n > 50. 
Since we have no theory that would determine the maxi- 
mum exponent n appearing in a superluminal dispersion 
relation, it would seem artificial at this stage to adjust 
n in order to achieve a sub-Planckian resolution of the 
horizon problem. 
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